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ABSTRACT 

We analyze publicly available void catalogs of the Baryon Oscillation Spectroscopic 
Survey Data Release 10 at redshifts 0.4 < z < 0.7. The first goal of this paper is to 
extend the Cosmic Microwave Background stacking analysis of previous spectroscopic 
void samples at z < 0.4. In addition, the DR10 void catalog provides the first chance 
to spectroscopically probe the volume of the Granett et al. (2008) supervoid catalog 
that constitutes the only set of voids which has shown a significant detection of a 
cross-correlation signal between void locations and average CMB chill. We found that 
the positions of voids identified in the spectroscopic DR10 galaxy catalog typically 
do not coincide with the locations of the Granett et al. supervoids in the overlapping 
volume, in spite of the presence of large underdense regions of high void-density in 
DR10. This failure to locate the same structures with spectroscopic redshifts may arise 
due to systematic differences in the properties of voids detected in photometric and 
spectroscopic samples. In the stacking measurement, we first find a AT = —11.5 ± 
3.7 fiK imprint for 35 of the 50 Granett et al. supervoids available in the DR10 volume. 
For the DR10 void catalog, lacking a prior on the number of voids to be considered in 
the stacking analysis, we find that the correlation measurement is fully consistent with 
no correlation. However, the measurement peaks with amplitude AT = —9.8 ±4.8 fiK 
for the a posteriori-selected 44 largest voids of size R > 65 h _1 Mpc that does match in 
terms of amplitude and number of structures the Granett et al. observation, although 
at different void positions. 

Key words: surveys - cosmology: observations - large-scale structure of Universe - 
cosmic background radiation 


1 INTRODUCTION 

Large-scale structures at low redshift leave their mark on 
the Cosmic Microwave Background (CMB) radiation pro¬ 
viding direct probes of the late ti me cosmic acceleratio n 
and the physics of Dark Energy (lAghanim et al] 120081) . 
In particular, large voids and clusters can imprint them¬ 
selves to the primary fluctuations of the CMB via phys- 
ical mechanisms call ed the Integrated Sachs-Wolfe effect 
(I Sachs fc Wolfel 1 19671. ISW) in the linear regime, and the 
Rees-Sciama effect ( Rees fe Sciamalfl968l . RS) on non-linear 
scales. The expected ISW correlation in the A CDM model 
is on the order of 0.1 41 K < |AT| < 1 yiK for typical voids 
dCai et al. IH), extending up to |AT| r; 20 ^tK for the 
largest observable superstruc t ures which are a l so the rarest 
(see e.g. ISzapudi et al.l d2014l ); lNadathur et alJ d2014l )). The 
contribution of the non-linear RS effects remains typically 
at the ~ 10% level compared to the linear expectation 
dCai et alj|2010h . However, the ISW and RS effects and their 


relative strength may be different in alternative cosmologies 
dCai et alj|2014h . 

The typical ISW and RS imprints are thus small enough 
to be immeasurable in practice. The traditional approach for 
detecting the weak signal is the angular cross-correlation 
measurement between galaxy density maps and the CMB. 
This detection strategy has been follo wed b y a series of 
studie s finding both margina lly (see e.g. iFrancis fc Peacock 
d20ld) . iKovacs et al. d2013h I and moderate ly sign ifican t 
(see e.g. m et al.l d2008l); iGiannantonio et al.l (12008 . 120121 ) ; 
IPlanck Collaboration et aid ( 2014al ) and references therein) 
ISW-like signals. Another approach is focussed on the largest 
structures in the density field, where t he ISW- RS effect is ex - 
pected to be the strongest. Foremost, I Granett et al.l d2008l ) 
created a catalog of supervoids and superclusterB (Gr08, 


http://ifa.hawaii.edu/cosmowave/supervoids/ 
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2 Andras Kovacs and Benjamin R. Granett 


hereafter) using the SDSS Data R e lease 4 (D R4) Mega-z 
photometric LRG catalog (IColIister et al.l[20071 ') with some 
additional area from DR6. They used the ZOBOV algorithm^ 
(lNevrinckll2008l 'l to identify the most prominent extremes of 
the large-scale density field. The superstructures were then 
used for stacking the CMB temperature centered on these 
locations, and to measure an average effect through a com¬ 
pensated top-hat filter. 

This rather simple statistic averages the CMB temper¬ 
ature centered on voids within a circular aperture r < R, 
and then subtracts the background temperature averaged 
in an equal-area con centric annulus with R < r < \plR. 

I Granett et alJ (l2008l l found a |AT| = 9.6 ± 2.2 /iK signal 
for their 100 most significant (> 3cr) superstructures us¬ 
ing an aperture size of R = 4°. This signal appears to be 
in 2cr tension with AC DM predictions, as pointed out 
in several follow-up studies using theory and simulations 


( Papai et all 20fd|_ 

20121; iFlender et al. 


Papai fe Szapudi 2010; Nadathur et al.l 


120131: iHernandcz-Monteagudo fe Smithl 


2013: lHotchkiss et al. 2015l ~ Aiola et al. 2015l b The expecta- 

tion for the stacked ISW signal remains at the |AT| ^ 2 fiK 
level for ~ 50 superstructures similar to Gr08 objects. 

Also, numerous additional tests have been performed 
to unc over possible s y stematic problems and statistical 
biases (|llic et all 120131 : IPlanck Collaboration et al] l2014al : 
ICai et all 20141 '). It was found that varying the number of 
the objects in the stacking, or using different filter sizes typ¬ 
ically lowers the overall significance. Otherwise the original 
Gr08 signal has survived every revision and remains a puz¬ 
zle. _ _ 

Additi onally , IPlanck Collaboration et al.l (l2014al f , 
ICai et aTI (l2014h . and |lh4_et_al] (|2013|) repea ted the CMB 
stacking analysis of Granett et al.l l l2008lj using com¬ 
plementary void catalogs ( Sutter et al. 201211 based on 
spectroscopic measurements at z < 0.4. These studies, 
however, did not report a highly significant detection of the 
ISW-like effect foun d in Gr08, e xcept the ~ 2cr evidence for 
a correlation in lCai et al. IiUHi). 

Due to the overwhelming contribution from cosmic vari¬ 
ance, ISW stacking measurements can be pron e to misinter - 
pretation due to the ‘look elsewhere’ effect llPeiri IHH). 
The measurement parameter space has a high dimension¬ 
ality when counting the choices made in the void catalogue 
selection and methodology details such as filter size, and the 
inherently weak s ignal can lead to the over interpretation of 
statistical flukes (iHernandez-Monteaeudo fe Smithll2013l l. 

In this work we do not formally carry out a blind analy¬ 
sis to mitigate the effects; however, we fix the measurement 
methodology using parameters determined externally. We 
test the robustness of the signal by varying these parame¬ 
ters, expecting that a true signal will be robust to perturba¬ 
tions in the catalogue properties or methodology. We com¬ 
pare the Gr08 void catalog to the Baryon Oscillation Spec¬ 
troscopic Survey (B OSS) DRIP C MASS and LOW-Z void 
catalogs provided by ISutter et al.l (|2014h . On one hand, we 
are able to probe a significant fraction of the Gr08 volume 
(surveyed with photo-z data) now with voids detected using 
spectroscopic redshifts. On the other hand, we extend pre¬ 
vious low-z DR7 void stacking measurements to the range 


2 http://skysrv.pha.j hu.edu/neyrinck/voboz/ 


of 0.4 < 2 < 0.7. Two distinct con clusions are possible: our 
stacking analysis could confirm the lGranett et aD (l2008l l de¬ 
tection for the first time with an independent void catalog, 
or the ISW(-like) signal could disappear. In any case, some 
puzzle will certainly remain for future analyses with BOSS 
DR12 and Planck DR2. 

The paper is organized as follows. Data sets, algorithms, 
and our observational results are presented in Section 2; the 
final section contains a summary, discussion and interpreta¬ 
tion of our results. 


2 DATA SETS AND MEASUREMENTS 
2.1 CMB data 


On the CMB side we use Planck’s SMICA0 map 
(IPlanck Collaboration et alJ l2014bl l with reso lution down¬ 
grade d to N s id e = 512 with the HEALPix llGorski et al.l 
I2005T) pixelization. We mask out potentially contaminated 
CMB pixels using t he WMAP 9-year ex tended tempera¬ 
ture analysis maslfl (iHinshaw et al.|[2013h at N s id e = 512 
to avoid repixelization effects of the N s ide = 2048 CMB 
ma sks provided b y Plan ck. It has already been p ointed out 
by iGran ett et all (|2008l) , and later confirm ed by [llic_et - ^il 


(20R), Planck Collaboration et ahl (l2014al l , and ICai et al 


( 20141) that the ISW-like cross-correlation signal detected at 
void locations is independent of the CMB data set when 
looking at WMAP Q, V, W, or Planck SMICA tempera¬ 
ture maps. We thus limit our analysis to the latest Planck 
SMICA sky map. 


2.2 Catalogs of cosmic voids 

We use public void catalogs by ISutter et al] (120141 ') where 
the authors ide ntified voids in BOSS DR10 spectroscopic 
galaxy samples (lAhn et al. 2014). The voids were identified 
with the ZOBOV too l llNeyrinckl2008ll within the VIDI0 frame¬ 
work (ISutter et al.ll20151 ). The void-finder ZOBOV is based on 
the watershed algorithm which builds a hierarchy of under¬ 
densities. We first remove voids with size R < 40 h _1 Mpc 
to cut the hierarchy. 

Supervoids encompassing smaller voids, walls, and fil¬ 
aments may be represented by multiple voids in a dense 
tracer catalog, or by one large void in sparsely sampled 
data or in the presenc e of photometric redshift uncertain¬ 
ties (IS utter et al.ll2014h . The galaxy number density for the 
CMASS catalog peaks at z=0.5 with mean inter-galaxy sep¬ 
aration L ~ 16 h _1 Mpc which rises to ~ 25 h _1 Mpc at 
z = 0.65, thus a lower cut of roughly twice this character¬ 
istic scale is a safe and reasonable choice to prune spurious 
void detec tions that would contribu te only noise to the mea¬ 
surement (ISutter et al.ll20fl l2015ll . 

Furthermore, these small a nd potenti ally spurious voids 
may occupy over dense regions. ICai et al.l (120141 ') tested this 
effect in mock catalogs. They suggest a lower size cut of 
R > 65 h^Mpc for a complete removal of potentially 
spurious voids at z me d « 0.43 in their Irgbright sample 


3 http://™.cosmos.esa.int/web/planck 

4 http://lambda.gsfc.nasa.gov/product/map/dr5/ 

5 http://www.cosmicvoids.net 
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Cold imprint of supervoids with Planck and BOSS DR10 3 


Denisty of DR10 voids vs. Gr08 void positions 



Figure 1 . Void positions of the Gr08 sample (gold) vs. DR10 void catalog (purple). Red points show a subclass of large DR10 voids of 
size R v > 60 h — 1 Mpc. The underlying = 32 HEALPix map is the DR10 void density sample smoothed with a a = 4° Gaussian. Due 

to the hierarchical organisation of the void catalogue, an underdense region may be split into a number of voids of various sizes in the 
catalogue. Dark blue colors indicate higher void abundance thus lower average projected density. The fact that the two void catalogues 
trace different structures points to systematic differences in the galaxy catalogues and void finding algorithms. 


(L w 38 h _1 Mpc), i.e. the DRY subsample that is possi¬ 
bly the most similar to the better sampled DR10 data we 
consider here. Although it is not possible to adopt such cuts 
for DR10 voids without proper simulations of the source 
catalogs, a rough comparison of the source densities of DR7 
and DR10 catalogs indicates that R. > 40 h _1 Mpc might be 
a reasonable cut. 

We also restrict our analysis to central voids to minimise 
possible contamina tions caused by the survey mask (see e.g. 
ISutter et al.l (l2014l l for details). 

These moderately conservative cuts remove ~ 65% of 
the voids from the DR10 catalog. Our approach is to probe 
the 0.4 < z < 0.7 redshift range thus we consider CMASS 
data in our analysis, with 13 extra voids from the LOW-Z 
4 sample at 0.4 < z < 0.45, i.e. a redshift range of the Gr08 
catalog not co vered by CMASS da ta. 

Following ISutter et al.l (120141 ') . we divide the result¬ 
ing DR10 void catalog into three subsamples; a combined 
CMASS 1 + LOW-Z 4 at 0.4 < 2 < 0.5 (56+13 voids), 
CMASS 2 at 0.5 < 2 < 0.6 (237 voids), and CMASS 3 at 
0.6 < z < 0.7 (172 voids). We analyze these catalogs both 
separately and jointly. We removed 6 voids from the anal¬ 
ysis, as their position was curiously outside of the rough 
Nside = 32 DR10 footprint by ~ 2°, indicating an inconsis¬ 
tency in the void catalog. We checked the effects of these 
objects on our main results and found no difference. 

6 We consulted Sutter et al. about this issue, and learned that 


We show our void sample together with the 50 Gr08 su¬ 
pervoids in Figure 1 (35 of the 50 Gr08 voids should be de¬ 
tectable in DR10, i.e. not masked out or residing close to the 
boundary). On average DR10 voids are smaller in angular 
and physical size than the Gr08 superstructures, due to the 
ability of resolving small-scale structures with spectroscopic 
redshifts. However, we found large regions of high density of 
DR10 voids, which typically do not overlap with the larger 
Gr08 supervoids. Thus the expected fragmentation of su¬ 
perstructures into smaller voids is not observed, but addi¬ 
tional large underdense regions appear in the spectroscopic 
data. This somewhat counter-intuitive finding means that 
a potential (and expected) ISW(-like) signal in the shared 
DR10-Gr08 volume is in this case carried by voids at distinct 
locations in the sky. 

The number densities of the two catalogues are of sim¬ 
ilar order: while the CMASS tracer inter-galaxy separation 
is fi -1 / 3 ~ 14—22 h _1 Mpc, the photometric LRG sample 
is ~ 2.2 times more dense with h -1 ^ 3 ~ 10 — 17 h _1 Mpc. 

The largest DR10 voids (minimum as large as the typi¬ 
cal Gr08 supervoids), however, are located closer to the Gr08 
supervoids on average. DR10 voids of size R > 60 h _1 Mpc 
are shown separately in Figure 1. Appropriate analysis of 
the properties of this subclass of voids is one of the main 
goals of this paper. 

these objects will be revised in a later version of their DR10 void 
catalog, and should be removed at the moment. 
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Figure 2. Filtered temperatures in re-scaled top-hats are shown as a function of their physical size. Crosses in the left panel mark 
LOW-Z 4 voids. Color bars indicate the redshifts of the voids, without any apparent trend or clustering in this parameter space. We note, 
however, that there is a slight extra average cooling for small CMASS 2 voids. Interestingly, CMASS 3 voids behave inversely showing 
hotter temperature differences on average for the smallest voids. Shaded regions mark out 2a statistical uncertainties scaled with the 
number of objects considered, while solid lines indicate the stacked temperature for a given subsample. 


2.3 Methods & Results 


Forem ost, we repeated the stacking analysis of lGranett et al.l 
(l2008h for the 35 supervoids available in the DR10 volume 
with constant R = 4° filter radius. The original sig na l of 
AT = —11.3 ± 3.1 /xK signal as measured bv lGranett et all 
(l2008lf now limited to the DR10 area is changed to AT = 
—11.5 ± 3.7 fj, K. We then expected to detect a similar sig¬ 
nal in the same physical volume with voids identified using 
spectroscopic redshift from the DR10 CMASS catalog. 

In ou r meth ods, we closely follow lllic et all (120131 ') and 
ICai et ahl t)2014h . We first measure average temperatures in 
the SMICA map at void locations using the comp e nsate d 
top-hat filter technique applied by iGranett et al.l (l2008h . 
W e further s cale t he filter b y angu lar si ze as advanced 
by Ilic et all ll2013l l , ICai et ahl (l2014| j , and iHotchkiss et al.l 
(l2015h . The same authors empirically found in data and in 
simulations that the optimal filter size to maximize S/N is 
~ 60% of the re-scaled void radius. The physical motivation 
behind such a scaling is the coincidence with the zero cross¬ 
ings of the void density profile and th e cumulative ISW -RS 
signal found in N-body simulations (ICai et al.l l2014h . We 
adopt this refinement in order to maximize the expected 
signal-to-noise in our tests. The resulting typical re-scaled 
filter radius is r mea n ~ 1-2° for all sub-samples. 

Our findings are presented in Figures 2 and 3, showing 
the void temperature measured as a function of void ra¬ 
dius and redshift. We compare the redshift distributions of 
DR10 and Gr08 voids in Figure 3. These plots show the typ¬ 
ical behavior of such top-hat filtered temperatures, as they 
contain large fluctuations for individual objects of both pos¬ 
itive and negative signs. However, there is no obvious excess 
clustering or other oddity in these parameter spaces. Two 
exceptions are the slight average shift to the negative side 
for R < 50 h^ 1 Mpc CMASS 2 objects, and the counterac¬ 
tive change at R < 50 h _1 Mpc for CMASS 3. Note that 
these voids should carry the lowest ISW-RS signal among 
the catalog, and their robustn ess is questiona ble due to the 
occurrence of voids in clouds (ICai et al.ll2014f ). 

We then average the filtered temperatures for the 478 
DR10 voids, sorted by radius. We estimated statistical un¬ 
certainties by analyzing 1000 Gaussian CMB simulations 
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Figure 3. Filtered temperatures in re-scaled top-hats are shown 
as a function of their redshift, covering the full range of Gr08 
voids. No meaningful trend is observable, as all sub-samples show 
similar distributions. The grey shaded region marks the 1 <j fluc¬ 
tuation ctat ~ 32 j/K for a single object, as measured using CMB 
simulations. 


generated with the HEALPix (iGorski et alj l2005l f synfast 
routin e using the Planck DR1 b est fit CMB power spec¬ 
trum (IPlanck Collaboration et al.l l2014bll. Gaussian simula¬ 
tions without considering instrument noise suffice because 
the CMB signal is domin ated by cosmic varianc e on the 
scales we consider (See e.g. IHotchkiss et ahl J2015I ')'). 

The ISW-RS signal expected in ACDM is so small 
that it is dominated by the uncertainty of the primary 
anisotropies even with stacking applied. More precisely, 
S/N < 0.4 for filters R < 2° was estimated bv lFlender et al.l 
(I2013T) . iNadathur et al.l ll2012ll pointed out that N sta ck ~ 
3000 supervoids could provide S/N ~ 2.5 for AT ^ —2 /xK, 
i.e. a signal that can be produced by the most extreme super- 
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Figure 4. Stacked CMB temperatures as a function of the num¬ 
ber of the objects considered. A physically motivated ordering of 
the voids by radius is applied, as the largest voids should leave 
the coldest imprints on the CMB. Dotted lines represent our re¬ 
sults when an additional weighting by void probability is applied. 
Orange dashes in panel 2 mark the errors obtained by measur¬ 
ing standard deviations for all filter sizes in simulations inde¬ 
pendently. In the top panel, we compare our results to the best 
signal-to-noise achieved by non-optimal rescaling strategy. 


structures in AC DM. Unfortunately, current observational 
capabilities cannot provide such a numerous catalog of voids. 

We adap t the error analysis for the stacki n g pre ¬ 
sent ed e.g. in IGranett et"ahl (| 20081 ). iFlender et alj (|2013l) . 
and ICai et all ( 20141) , and compare two error estimators. 
First, we repeated the stacking analysis 1000 times varying 
the CMB realization and fixing the position and scaling of 
the 56 top-hat filters from the CMASS 1 data. We com¬ 
pare this against a simpler approach in which the variance 
is computed for a single filter (assumed to be independent) 
and rescaled by the number of voids in the catalogue. The 
mean variance is estimated from 1000 realizations but on 
each realization the filter size is randomly drawn from the 
size distribution. This gives a mean standard deviation of 
(tat ~ 32 fiK and we find that this error does not depend 
on the redshift bin. We then scale the variance by the num¬ 
ber of objects in the stack: o^T{N sta ck) = 32 fiK/y/N/tack- 
This approach simplifies the computation and tests revealed 
that the two estimators agree at the per cent level based on 
the standard deviation measurement. The CMASS 1 panel 
of Figure 4 contains a comparison of the error bars obtained 



Figure 5. The probability of finding at least one spurious signal 
at the negative extreme with S/N > 2 is plotted as a function of 
N st ac k in the Gaussian noise model. The inset shows the prob¬ 
ability distribution of maximum S/N for the case of all voids 
considered in the stacking. 


with the two methods showin g good agreem e nt. Th is finding 
is in agreement with that of lGranett et al.l (120081 1 confirm¬ 
ing that the measurements with t he compensated filter may 
be considered to be independent. IGranett et al.l (l2008l l also 
pointed out that errors obtained by drawing random points 
to a given CMB map, and errors measured with fixed void 
positions varying the CMB realization agree at the ~ 2% 
level. 

The mean error we obtained is higher than the 
(tat ~ 22 uncer tainty for a single supervoid found by 

IGranett et al.l (l2008l l. due to larger CMB fluctuations at 
smaller scales. Note that this higher noise level prevents high 
signal-to-noise measurements using relatively small voids. 

This procedure allows us to estimate the significance of 
a measurement given a fixed number of voids Nstack- How¬ 
ever, with uncertainties in the properties of the voids and 
the origin of the signal we do not have a strong prior on the 
number to average in the analysis. Using only the largest 
voids we do not have the precision to measure a weak sig¬ 
nal, while adding the smaller voids we can wash-out a signal 
if it exists._ 

To address this situation, IGranett et al.l (120081 ) define a 
cutoff based upon void probability , taking all voids with de¬ 
tection signficance > 3 a. Further, ICai et akl (|2014h applied 
a weighting as a function of void probability; however, they 
found no significant difference in the measured signal. We 
also test this weighting scheme here, as shown in Fig. 4 and 
discussed below. 

Without using a prior as to which voids to consider 
in the measurement we are led to examine the signal as a 
function of N s t ac k ordered by void size, as in Fig. [4] The 
temperature measured after averaging N sta ck structures ap¬ 
pears as a random walk and, given the null hypothesis of no 
correlation, we may estimate the probability that the tra- 
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6 Andras Kovacs and Benjamin R. Granett 


jectory crosses a given threshold for detection. We compute 
this probability in Monte Carlo fashion assuming the Gaus¬ 
sian noise model. We generate a vector of Ntot Gaussian 
distributed values with a = 32 /iK. We then check if the cu¬ 
mulative sum crosses a given significance threshold at any 
point. This is repeated and we keep count of the fraction 
of runs that give a significance above 2cr. The result after 
100000 runs is shown as a function of N stac k in Fig. [5] We 
find that given a catalogue of 478 voids, the odds of finding a 
2 a signal in the cumulative temperature measurement given 
the null hypothesis is 24%. 

We now consider the stacked temperature measure¬ 
ments shown in Fig. 4. 

(i) Considering the combined sample (Fig. 4., top panel), 
the signal fluctuates around the la level. A cold imprint 
peaks with an amplitude of —9.8 ± 4.8 /iK or ~ 2<r for the 
44 largest voids with sizes R > 65 h -1 Mpc. However, this is 
unremarkable given that the probability of finding such an 
extreme somewhere in the cumulative stacked temperature 
with a total of 478 structures is ~ 24%. 

(ii) However, the Gr08 sample contains 50 significant su¬ 
pervoids, thus in terms of the largest ^50 fluctuations the 
two catalogs seem to agree with each other even if the struc¬ 
tures differ in size and are not at the same positions. 

(iii) Counting all 478 voids down to R = 40 h _1 Mpc the 
signal sharply approaches to zero and becomes AT « 1.2 ± 
1.5 /iK. Note that the overall measurement, including the 
peak, is unchanged when weighting by void probability. 

(iv) Applying a presumably non-optimal filter scaling of 
40% of the void radius, we found —5.8 ± 3.0 /rK or ~ 1.9<r 
for the 118 largest voids with sizes R > 55 h -1 Mpc. This 
signal, shown by the error bar in the top panel of Figure 4, is 
the largest measured using different rescaled filter sizes. Note 
that this signal is also observed at relatively large void radii, 
and it is consistent with our measurement at Nstack = 118 
using R/R v = 0.6. 

(v) The combined CMASS 1+LOW-Z 4 sample shows a 
signal fluctuating inside the la level, resulting in a final value 
of AT Ri 1.3 ± 3.9 /iK for the full sub-catalog with 69 void 
members. Note the effect of the probability weighting which 
reduces the amplitude. 

(vi) The CMASS 2 sample with the largest number of 
voids contributes most strongly to the combined sample, 
thus the signal is similar to that described above. After fluc¬ 
tuating inside the ler expectation and reaching the ~ 1.5<r 
level or AT « —10 /iK at R > 65 h~ 1 Mpc, the overall signal- 
to-noise with all voids included is AT « —1.7 ± 2.1 /iK for 
the 237 CMASS 2 objects. The weighting by void probability 
results in a more effective convergence to zero. 

(vii) The CMASS 3 temperature signal fluctuates around 
the la level. Adding smaller scale voids the signal becomes 
positive and results in a curious positive AT « 5.3±2.5 /iK. 
However, the galaxy number density is lowest in the CMASS 
3 bin, so a larger fraction of the smallest voids may be spuri¬ 
ous compared with the other more densely sampled redshift 
bins. This intuition is verified by the probability weighting 
test, that lowers the significance of the spurious positive sig¬ 
nal, resulting in AT w 2.0 ± 2.5 /iK. 




^stavk 


Figure 6. Stacked CMB temperatures for voids of size R v > 
60 h — -'Mpc ordered by radius. In the bottom panel, we show the 
signal-to-noise ratio of our measurement for different re-scaling 
parameters. In the top panel, we compare various measurement 
strategies for case R/R v = 0.6, i.e. the best filter size obtained 
in both simulation and data. The maximal signal-to-noise is ob¬ 
served at the lower size limit of R v > 65 h —1 Mpc. The signal we 
detected is comparable to the Gr08 observation (shown in blue), 
although less significant. 

2.4 Imprint of DR10 supervoids 

We have an indication that the largest structures in DR10 
may produce an imprint on the CMB characterised by the 
—9.8 ± 4.8 /iK or ~ 2cr signal that we measured for voids of 
size R > 65 h” 1 Mpc. Although on its own the occurrence of 
the signal is unremarkable in the light of our Monte Carlo 
tests, it does match the Gr08 measurement in terms of am¬ 
plitude and number of structures. Thus this result warrants 
further inspection. Our findings, based on additional sys¬ 
tematic tests are shown in Figure 6, and summarized below 

(i) the lower panel of Figure 6 demonstrates that the 
R/R v = 0.6 scaling gives the highest S/N, as found in sim¬ 
ulations and DR7 data. 

(ii) the shape of the dotted purple curve in the upper 
panel indicates that this signal is rob ust against wei ghting 
by void probability, in agreement with lCai et al.1 (|2014l ). 

(iii) we randomly changed the void radius by +10% or 
— 10% for all voids, and the signal nearly disappeared, as 
shown in blue. 

(iv) we artificially removed large-scale fluctuations (i.e. 
much larger than the void size) from the SMICA CMB map 
using HEALPix routines, and repeated the stacking. Without 
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Figure 7. A comparison of filtered CMB temperatures as a func¬ 
tion of the physical size of voids. Shaded regions mark 2cr fluc¬ 
tuations, while solid lines show the measurements by Gr08 and a 
result obtained using our final sub-catalog. 


(SZ) effect have already been excluded by Gr08 with conser¬ 
vative masking and probes of CMB color dependence further 
const rained by Planck analyses (iPlanck Collaboration et ali 
l2015lh Another uncertainty in the superstructure catalog, in¬ 
troduced by a z ~ 0.05 photometric redshift errors of Mega-z 
LRGs, have not been investigated. 

Figure 7 shows a comparison of filtered temperatures 
for Gr08 and DR10 voids as a function of their radius. In 
summary, the main properties are as follows. 

(i) larger fluctuations are observed in DR10 results due 
to smaller filter size. 

(ii) proportionally more R > 65 h -1 Mpc (super)voids 
were found in Gr08 than in DR10, possibly due to systematic 
differences in void detection and source catalogs. 

(iii) in Gr08, the largest supervoids leave an ever stronger 
stacked imprint in the CMB. 

Large underdense regions of diameter ~ 130 h _1 Mpc, 
perhaps better thought of as long wavelength fluctuations, 
can contain smaller voids, walls, and filaments. Therefore, 
to identify these supervoids in galaxy surveys it will be nec¬ 
essary to modify void detection algorithms. The voids ef¬ 
fectively resolved by spectroscopic redshift surveys may be 
grouped together and unified to identify under-densities on 
the largest scales. Focusing on the class of supervoids would 
further simplify the interpretation of results by limiting the 
selection effects arising from a posteriori choices and illu¬ 
minate the possible connection between the apparent ISW 
excess and the density field on the largest scales. 


£ < 10 modes, the map has changed, and errors decreased by 
~ 3%, but the filtered temperatures remain almost the same 
(shown with green dashes ). This fin d ing is again consistent 
with the DR7 analysis bv ICai et al.l (120141 ') . 

(v) despite the different stacking methods and void prop¬ 
erties, the signal we found is in lcr agreement with Gr08, 
although less significant due to larger CMB fluctuations at 
smaller angular scales (Gr08: constant R = 4° filter, DR10: 
R/R v =0.6 re-scaling) 

Therefore, we see that the our measurement appears 
to be robust against physically motivated changes in the 
analysis, while the weak signal disappears when non-optimal 
techniques are applied. The exception is when the void radii 
were perturbed by 10% and we find a reduction in the signal. 
An inherent uncertainty is expected in the void radii due to 
sampling error and variation in void shape, so this may point 
to a fragility of the signal. 

Keeping in mind that the ISW-RS signal expected in 
A CDM remains below AT ss —2.0 pK even for the 50 
largest posteriori-selected supervoids we also lack the the¬ 
oretical motivation for this signal. 

While the 4.4cr signal for the posteriori-selected, 50 most 
significant Gr08 superstructures appears to be too large for a 
plausi ble primordial fluctuation. IPlanck Collaboration et al.l 
(1201511 explicitly pointed out that the lack of a strong T ~ E 
correlation in the CMB polarization data by Planck in the 
location of the Gr08 structures provides evidence that the 
effect is not a primordial fluctuation. Also, systematic ef¬ 
fects are obvious suspects, but explanations such as spurious 
correlations caused by stellar contamination, extragalactic 
radio sources, or contamination via the Sunyaev-Zeldovich 


3 CONCLUSIONS 


We probed the volume of the iGranett et al.l (|2008l l super- 
vo id catalog with the BOSS DR10 void catalog provided 
by lSutter et al.l (120141 '). Our principal aim was to revisit the 
strong ISW(-like) signal found in Gr08 with a catalog that 
probes the same density f ield. We p r uned the DR10 catalog 
followi ng the protoc o l of Cai et al.l (120141 4 and the sugges¬ 
tion of ISutter et al.l (|2014i l for removing the smallest and 
least reliable voids which are also expected to produce the 
smallest ISW-RS signals in A CDM. The voids identified 
with spectroscopic redshifts in DR10 are smaller than the 
Gr08 structures traced with photometric redshifts. Even so, 
we find that the Gr08 supervoid positions do not coincide 
with regions abundant in DR10 voids which indicates that 
the void finders are sensitive to different structures. This 
situation merits further study to understand the systematic 
differences between voids identified in spectroscopic versus 
photometric samples. 

We measured the stacked CMB temperatures for the 
35 Gr08 supervoids in the DR10 footprint using the orig¬ 
inal filter size. The AT = —11.5 ± 3.7 pK signal that we 
found is consistent with the Gr08 measurement. We then 
performed a stacking analysis using 478 DR10 CMASS and 
LOW-Z voids. In general, no significant imprint was detected 
in the presence of large cosmic variance. However, we found 
a AT = —9.8±4.8 fiK or 2cr signal for the largest 44 voids of 
size R > 65 h _1 Mpc in the combined sample. We note that 
this detection in itself is unremarkable without a strong prior 
on the number of sources or minimum void size used in the 
analysis, as we found that the probability of finding a 2er sig- 
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nal somewhere in the cumulative temperature measurement 
with 478 voids is p ~ 24%. However, the measurement is of 
interest considering the detection based on the Gr08 sample. 
We examined how robust the measurement is to variations 
in the methodology. Varying the filter size, we found that 
the signal is maximised when using a filter radius found to 
be optimal in N-body simulations. 

Our results highlight that ISW detections with the 
stacking protocol strongly depend on the properties of the 
tracer population and the void finder. While the effect of 
photo-z errors on the performance of ZOBOV ha ve not been 
tested, it has been emphasized bv ISutter et al.1 (l2014lf that 
masking and the density of the tracer population strongly 
affects the resulting void catalogs. Further systematic un¬ 
certainties remain in how the hierarchy of voids is cut and 
the location of void centers. 
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